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ABSTRACT 

On the Sun, the rotation periods of individual sunspots not only trace the latitude 
dependence of the surface rotation rate, but also provide clues as to the amount of 
subsurface fluid shear. In this paper we present the first measurements of stellar dif- 
ferential rotation made by tracking the rotation of individual starspots with sizes 
comparable to the largest sunspots. To achieve this we re-analyse four sequences of 
densely-sampled, high signal-to-noise echelle spectra of AB Doradus spanning several 
stellar rotations in 1996 December. Using spectral subtraction, least-squares decon- 
volution and matched-filter analysis, we demonstrate that it is possible to measure 
directly the velocity amplitudes and rotation periods of large numbers of individual 
starspots at low to intermediate latitude. We derive values for the equatorial rota- 
tion rate and the magnitude of the surface d ifferential rotation, both of which a rc in 



excellent agreement with those obtained by Donati fc Collier Cameron (1997) from 
cross-correlation of Doppler images derived a year earlier in 1995 D ecember, and with 
a re-analysis of the 1996 data by the method of Donati ct al. (2000)} The differences be- 



tween the rotation rates of individual spots and the fitted differential rotation law are 
substantially greater than the observational errors. The smaller spots show a greater 
scatter about the mean relation than the larger ones, which suggests that buffeting 
by turbulent supergranular flows could be responsible. 

Key words: stars: activity - stars: imaging ~ stars: individual: AB Dor - stars: 
rotation - stars: spots 



1 INTRODUCTION 

Observational studies of differential rotation on rapidly ro- 
tating st ars have recently been published by a number of 



author s. Donati & Collier Cameron (1997) 



and 



Donati et al 



(1999) applied a latitude- by- latitude cross-correlation analy- 



sis to pairs of surface images of the rapidly rotating K dwarf 
AB Doradus secured several days apart. They found a solar- 
like rate of surface shear, with the equatorial region rotating 
faster than higher latitudes and a beat period of about 110 
days between the equatori al and high-latitude rotation pe- 
riods. Barnes et al. (2000} obtained a similar result for the 



young K dwarf PZ Tel. [Donati et al. (2000)| used a para- 



metric imaging method to establish the differential rotation 
rate o n the pre-main sequ ence G star RXJ1508-4423 (= LQ 
Lup). Petit et al. (2001) measured the differential rotation 
of the K subgiant prim ary of the RS CVn bina ry HR 1099 
by this same me thod. Vogt et al. (1999) and Strassmeier 
& Bartus (2000) examined the evolving starspot distribu- 
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tion on HR 1099 from long sequences of Doppler images, 
and found evidence of a gradual poleward drift of large spot 
groups on timescales of months to years. 

Studies based on Doppler images suffer from a number 
of drawbacks for this kind of work. When the sampling of 
the line profile around the stellar rotation cycle is sparse, 
the problem of mapping the stellar surface-brightness dis- 
tribution becomes an ill-conditioned inverse problem, neces- 
sitating a regularised fitting approach. The regularising 
function imposes a penalty on the growth of fine-scale struc- 
ture in the stellar image. This has the well-documented ef- 
fect of biassing the centroids of low-latitude features in the 
images away from the stellar equator, and may even sup- 
press them altogether. It it likely that differential rotation 
measurements based on these methods may suffer some dis- 
tortion as a result. 

Several recent Doppler imaging studies of rapidly ro- 
tating young main-sequence and pre-main sequence stars 
have, however, employed a much denser sampling strategy. 
With echelle spectrographs and large-format CCD detec- 
tors it is possible to monitor the profiles of thousands of 
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lines simultaneously in continuous sequences of exposures 
with a time resolution of a few minutes. Signal-enhancernent 
techniques such as least-squares deconvolution (LSD: |3o-| 
nati et al. 1997) stack up the profiles of thousands of pho- 
tospheric lines to give composite line profiles with signal-to- 
noise (S/N) ratios of several thousand. 

These densely-sampled series of composite line profiles 
delineate clearly the trajectories of line-profile distortions of 
many small spots. They bear a superficial resemblance to 
time-series spectra of rapidly rotating, early-type pulsating 
stars. Many of the analytic techniques that have been used to 
identify non-radial pulsation modes in these objects operate 
directly on the data, without recourse to large-scale inver 
sion techniques dKcnnclly, Walker fc Mcrryficld 1991 Wad ^ 



199i; B ghriivers fc Tcltiug 1991]|; Jaukov, ,Janot-Pacheco & 



Leister |2000|; [Tehing, Abbott fc Schrijvcrs 20011 ). Since the 
trajectories of starspot signatures in the trailed spectrogram 
of a rotating star are sinusoidal, it should be possible to use 
analogous forward-modelling techniques to identify individ- 
ual spots and determine their properties. 

The purpose of this paper is to apply matched-filter 
analysis to the problem of detecting and tracking individ- 
ual unresolved starspots in densely-sampled time sequences 
of high S:N line profiles. Section 2 presents a re-analysis 
of archival echelle spectra of AB Doradus observed in 1996 
December, from which Doppler images have previously been 
published. Section 3 develops a matched-filter methodology 
based on a similar technique used by Collier Cameron et 
al (1999, 2001) to search for faint reflected-light signatures 
of extra-solar planets. In section 4 this technique is used to 
measure the rotation periods of individual spot signatures 
present in each of these data sets, as a function of their 
rotational velocity amplitudes. In section 5 the results are 
compared with those predicted by differential rotation mod- 
els derived from previous Doppler imaging studies. 
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Figure 1. The upper panel shows a sample least-squares decon- 
volved profile of AB Dor, computed using a list of stellar photo- 
spheric lines only (upper), and augmented with a list of telluric 
lines (lower). The asymmetries and bumps in the profile are the 
signatures of starspots. The difference between these two profiles 
(lower panel) shows a clear peak near zero velocity, produced by 
telluric H2O and O2 lines. 



2 DATA PREPARATION 

The datasets employed here were secured using the 
SEMELPOL dual-fibre polarimeter feed to the UCL echelle 
spectrograph (UCLES) at the 4 .2-m Anglo- Austra lian Tele- 
scope on 1996 December 23-29 (Donati et al. 199£). Readers 
are referred to the above paper for the journal of observation 
and details of the instrument configuration, and the initial 
extraction of the echelle spectra from the raw data frames. 

The starspot signatures were isolated using a spectral 
subtraction and least-squares deconvolution procedure. The 
first task was to place all the spectra on a common velocity 
scale. This was achieved via th e least-squares deconvolution 
method of Donati et al. (1997) . The essence of this technique 



is that the observed spectrum is treated as the convolution of 
an "average" line profile and a list of delta functions at the 
wavelengths of the known photospheric lines derived from 
a Kurucz model. The lines are weighted according to their 
central depths, and the "average" profile is computed using 
the method of least squares, such that it gives an optimal 
fit to the spectrum when convolved with the line list. 

If the line list is augmented to include a list of known 
telluric lines in the observed wavelength range, the decon- 
volved profile develops a small, narrow pseudo-absorption 
bump near zero velocity. This feature becomes more appar- 



ent when we subtract the we subtract the deconvolved profile 
computed using the stellar line list only, as shown in Fig. ^ 
By fitting a gaussian to the central peak of this residual 
profile, we are able to determine the zero-point offset in- 
troduced by systematic errors in the wavelength calibration 
(using Th-Ar arc spectra taken at the start and end of each 
night) and by any instrumental flexure during the night. 
This offset is computed and applied to the zero-point of 
the velocity scale for all subsequent deconvolutions, thereby 
ensuring that the deconvolved spectra are accurately posi- 
tioned in the observer's rest frame. 

Telluric-line spectra derived from the rapidly-rotating 
B star HR 3084 were then scaled and subtracted from the 
individual echelle spectra. The individual spectra of AB 
Dor were then summed to create a time-averaged spectrum 
for each night of observation. This template spectrum was 
then scaled and subtracted from each individual spectrum 
in turn, using the spectral alig nment technique described by 
Collier Cameron et al. (2001)| . 

This left a residual spectrum consisting mainly of ran- 
dom noise, on which the residual line-profile distortions 
caused by stellar spots were superimposed. We deconvolved 
this residual spectrum with the stellar line list, this time ap- 
plying the heliocentric velocity correction computed for the 
mid-time of the exposure. This yielded profiles from which 
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Figure 2. Time-series spectra of deconvolved residual starspot signatures on tlic nigiits of 1996 December 23, 25, 27 and 29. Tiie 
starspot signatures appear as brigiit pseudo-emission trails crossi ng tlie stellar rota tion profile from negative to positive radial velocity. 
The rotation phases are computed according to the ephemeris of Innis et al. (1988), HJD = 2444296.575 -I- 0.51479i?. 



the average line profile had been removed, effectively isolat- 
ing the starspot signatures and placing them accurately on 
a heliocentric velocity scale. 

Any residual telluric features produced by misalign- 
ment of the telluric template and the target spectrum were 



removed at this stage using principal-component analysis 
(PC A) of the night's ensemble of spectra. Removal of the 
first two principal components was sufficient to eliminate 
the residual telluric lines without significantly affecting the 
starspot trails. A detailed description of the use of PCA 
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for removal of correlated noise features is given by [Gollier 
Camerc n et al. (2001). In practice this only proved necessary 



for the spectra observed on 1996 December 29. 

The resulting time-series of deconvolved starspot resid- 
uals are shown in Fig. ^ They show clearly the signatures 
of individual starspots, following sinusoidal paths about the 
stellar centre-of-mass velocity as they cross the visible hemi- 
sphere of the star. When the same range of rotation phases 
is observed on successive nights, the same features are seen 
to recur at more or less the same rotation phase. As Donati 
et al (1997, 1999) have shown, however, spots at high lati- 
tudes have slightly longer rotation periods than spots near 
the equator. High-latitude spots also have smaller velocity 
amplitudes, since they lie closer to the stellar rotation axis. 
Since we are interested in measuring the rotation periods 
and latitudes of individual spots in order to delineate the 
differential rotation pattern, we need to develop an optimal 
method for determining these two quantities directly from 
the trailed spectrograms. 



3 MATCHED-FILTER ANALYSIS 

The spectral signature of an unresolved spot on the stellar 
surface is closely approximated by a gaussian having the 
same full width at half maximum intensity (FWHM) as the 
local specific intensity profile of the spectral line concerned. 

The gaussian profile varies in amplitude and radial ve- 
locity according to the expression: 



giv,t\(l),K,n) 



/(7) 



exp 



K sin i sin ip — vq 
Ave 



Here vo is the radial velocity of the star's centre of mass. 
A spot located at latitude 9 on a spherical star lies at a 
distance J?^ cos 6 from the stellar rotation axis. If the stellar 
axis is inclined at an angle i to the line of sight, the spot's 
radial velocity will vary sinusoidally about vo with angular 
frequency Q{9) and amplitude K, where 



K = n{e)R^, cose sin i. 



(1) 



At any time t, a spot at position (0, 6) will be displaced 
from the observer's meridian in the direction of rotation by 
an angle 

The stellar longitude </!> increases in the direction of rotation, 
and is measured from the observer's meridian on the stellar 
surface at an arbitrary time to. The fiducial epoch to is cho- 
sen to be near the middle of the sequence of observations, 
to eliminate spurious correlations between </> and il. 

While the spot is on the visible hemisphere its apparent 
area will be foreshortened by a factor 

cos 7 — cos i sin 9 + sin i cos 8 costp. 

Here 7 is the angular separation of the spot and the observer 
subtended at the centre of the star. 

The flux blocked by the spot is also subject to limb dark- 
ening, which is assumed to have a linear dependence on cos 7 
with a limb-dark ening coefficient u — 0.77. This value, fro m 



the tabulation of Diaz-Cordoves, Claret & Gimenez (1995) 



inverse-variance weighted mean wavelength A = 5475A of 
the spectral lines used in the deconvolution. 

The overall reduction in signal strength for a spot at 
foreshortening angle 7 relative to its strength when viewed 
directly from above (7 = 0) is 

fil) = cos 7(1 — It + MCOS7). 

The gaussian width parameter Ave is related to the line 
FWHM by 

FWHM 
Ave = 



The spectrum of the slowly rotating K dwarf Gl 176.3 has 
been used in several previous Doppler imaging studies of AB 
Dor as a spectral template whose line proflles are represen- 
tative of the local line profiles on AB Dor in the absence 
of rotational broadening. A gaussian fit to the least-squares 
deconvolved profile of Gl 176.3 yielded Ave = 5.53 km s~^. 

The matched filter g{v,t\(j}, K,Q,) can be computed on 
a two-dimensional map with the same dimensions as the ob- 
served time-series spectrum, and fitted to the data via an 
optimal scaling. The scale factor W by which g is multiplied 
to give the best fit to an individual spot trail, is the equiv- 
alent width of the bump in the line profile that would be 
produced by the spot if it were located at the centre of the 
disc. It is proportional to the area of the spot, and related to 
the intensity ratio of the local continuum in the photosphere 
and the spot, the ratio of the line equivalent widths in the 
photosphere and the spot, the overall spot coverage, and the 
limb-darkening coefficient. 

3.1 Second-derivative profiles 

The matched-filter fitting procedure described above suffers 
from the disadvantage that although the residual starspot 
signatures are approximately gaussian, they appear super- 
imposed on a sloping, non-zero background that changes 
slowly from one spectrum to the next (Fig. ^). This means 
that W cannot be determined using a simple gaussian fit. 

Instead, we compute a numerical second derivative to 
pick out the "ridge line" traced out by a spot in the trailed 
spectrogram. The second derivative of each spectrum s in 
the observed time-series is defined as: 



-2, J 



+ 2,3- 



The choice of a two-pixel offset is a compromise between 
computing the derivatives over a small velocity range to 
avoid loss of resolution, and maximising the baseline over 
which the derivatives are computed to reduce noise. The 
smallest resolvable separation between neighbouring spot 
features appearing in a single line profile is roughly 5 pixels 
(Fig. which suggests that a two- pixel offset is optimal in 
this respect. 

We treat the matched filter in the same way: 



2/i,i + .fi+2,j 



is appropriate for an early-K photosphere at the depth- and 



Here fi^j is the contribution of the gaussian matched filter 
to the ith pixel of the jth profile in the time-series, which 
spans the velocity range from Wi_i/2 to Wi+1/2: 

fi,i = / 9ii'i,tj\(t),K,n)dv. 

•^■",-1/2 
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Figure 4. The right-hand subpancls show a function of rotational velocity amplitude K and phase, for four values of the rotation 

period: (a) 0.5125 d; (b) 0.5135 d; (c) 0.5145 d; (d) 0.5155 d. The corresponding trailed spectrograms in the left-hand subpanels show all 
four nights' data merged and phased using these rotation periods. A dark, oval peak in the right-hand panel corresponds to a bright 
starspot trail (white) crossing the centre of the line profile in the left subpanel at the same rotation phase. The velocity amplitude K in 
the right-hand panel is proportional to the radial a<;celeration of the spot signature as it crossed the centre of the profile. The definition 
of fine-scale features in the trailed spectrum is very sensitive to the period on which the data are folded. 



The scale factor W is then computed from these second- 
derivative maps via the optimal scaling: 

^_ Ef,j a:i,jP»,j/Var(gj,j) 
Ei.,PL/Var(xi,,) • 

Here crf j is the variance associated with the data point Sij, 
so 

Var(a;-ij) = cr^-2,J + 'ialj + (j1+2,j- 

Note that the scale factor W is the equivalent width in ve- 
locity units of the starspot bump that would be seen in the 

deconvolved profile if the spot wore observed at the centre 
of the stellar disc, and so is expressed in km s^^. 



4 RESULTS 



4.1 Spot parameter estimation 

The optimal fitting procedure determines both W and the 
associated badness-of-fit statistic 

2 _ {Xi,j - Wpi,jf 

for a given set of matched-filter parameters {(f>,K,Q). To 

identify spot signatures and measure their parameters, we 
construct a sequence of maps on a (0, K) grid, for a range 
of rotation periods P = 2tt/Q. 

The sequence of maps is stacked to form a data cube, 
W{(j), K.Q.) or Ax'^{4>, K,Q,) which is then searched for lo- 
cal maxima in all three dimensions. Inevitably some of the 
weaker local maxima are spurious noise features, so it is nec- 
essary to impose a cutoff on W (or better, on Ax^) below 
which candidate spots are rejected as being probable noise 
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Figure 3. Sample deconvolved profiles from the trailed spectro- 
gram for 1996 December 25 between phases 0.6 and 0.7 show 
closely-spaced peaks caused by marginally-resolved spots, sepa- 
rated by as little as 5 pixels. The residual intensity units are given 
as a fraction of the mean continuum intensity. 



Table 1 . Parameters derived from matched-filter analysis for the 
24 candidate starspot signatures in phase range 0.50 to 1.00 on 
1996 December 23, 25, 27 and 29. The six asterisked features are 
aliases. 



Phase 



0.5099 ± 
0.5471 ± 
0.5504 ± 
0.5514 ± 
0.5788 ± 
0.5878 ± 
0.6258 ± 
0.6284 ± 
0.6401 ± 
0.6737 ± 
0.6919 ± 
0.7252 ± 
0.7594 ± 
0.7878 ± 
0.8129 ± 
0.8134 ± 
0.8247 ± 
0.8788 ± 
0.9139 ± 
0.9327 ± 
0.9365 ± 
0.9580 ± 
0.9580 ± 
0.9687 ± 



0.0001 
0.0002 
0.0007 
0.0003 
0.0001 
0.0004 
0.0015 
0.0003 
0.0003 
0.0001 
0.0005 
0.0003 
0.0002 
0.0002 
0.0006 
0.0004 
0.0007 
0.0002 
0.0002 
0.0003 
0.0003 
0.0004 
0.0007 
0.0003 



K 

(km s"-*-) 



73.29 
87.32 
59.03 
83.41 
83.47 
86.79 
38.10 
80.67 
82.32 
81.48 
64.77 
81.30 
86.06 
81.04 
67.23 
91.01 
51.14 
76.77 
89.35 
85.78 
70.46 
81.67 
60.90 
85.76 



±0.07 
±0.34 
±0.60 
±0.29 
±0.04 
±0.19 
±0.47 
±0.18 
±0.09 
±0.12 
±0.36 
±0.17 
±0.10 
±0.09 
±0.66 
±0.16 
±0.14 
±0.14 
±0.11 
±0.46 
±0.25 
±0.12 
±0.52 
±0.06 



P 

(days) 



0.513964 
0.513299 
0.513966 
0.515664 
0.513597 
0.516632 
0.513507 
0.513437 
0.516597 
0.513839 
0.514308 
0.513334 
0.513738 
0.513662 
0.513837 
0.513032 
0.515061 
0.513847 
0.513543 
0.515426 
0.514310 
0.515077 
0.515384 
0.513512 



± 0.000023 
± 0.000037 
± 0.000088 
± 0.000036 
± 0.000020 
± 0.000063 
± 0.000371 
± 0.000055 
± 0.000052 
± 0.000019 
± 0.000083 
± 0.000047 
± 0.000037 
± 0.000026 
± 0.000065 
± 0.000066 
± 0.000092 
± 0.000026 
± 0.000035 
± 0.000057 
± 0.000037 
± 0.000055 
± 0.000081 
± 0.000033 



W 

(km s"-*-) 

0.00655 

0.00297 

0.00171 

0.00298* 

0.00596 

0.00176* 

0.00155* 

0.00267 

0.00196* 

0.00715 

0.00244 

0.00254 

0.00297 

0.00413 

0.00181 

0.00158 

0.00188 

0.00547 

0.00275 

0.00240* 

0.00369 

0.00201* 

0.00254 

0.00325 



features. An appropriate cutoff is found by computing the 
{(j), K) maps over a range of K values somewhat greater than 
the stellar wsin j, and tuning the cutoff to reject all features 
with K > V sin i. 

Note also from Fig. ^that good repeated phase coverage 
was obtained on at least 3 of the 4 nights only for features 
crossing the centre of the line profile in the phase range from 
0.50 to 1.00. The search was therefore restricted to this phase 
range. 

This yielded the list of candidate spots in Table |l|. The 
parameter values are listed together with their formal la 
errors, derived using 

and similarly for and P. The la errors on the spot ve- 
locity amplitude and period are consistent with the values 
expected from simple considerations: if the radial velocity 
of a spot bump can be determined to a precision of order 
0.1 pixel (300 m s~^), the uncertainty in the time at which 
it crosses the observer's meridian is of order 20s. Over a 6- 
day (12-rotation) baseline, this gives an uncertainty of about 
0.00003 day in the period determination. Determining the 
radial acceleration - and hence K - involves measuring the 
change in spot velocity over the course of roughly one-sixth 
of the rotation cycle, which again involves velocity and tim- 
ing measurements that are uncertain at the 0.1 km s~^ and 
20s levels respectively. The resulting uncertainty in K ought 
therefore to be of order 0.1 to 0.2 km s^^ for a well-observed 
spot. 

We explored the effects of possible errors in the limb- 
darkening coefficient by repeating the matched-filter anal- 



ysis using a more solar-like u = 0.6 rather than u — 0.77. 
Given that the time-averaged spectrum is subtracted prior 
to deconvolution, the principal effect of decreasing the limb- 
darkening coefficient is to decrease the contrast between the 
strength of a spot signature viewed at the centre of the disc 
and near the stellar limb. As expected, we found that the 
lower value of u yielded a slightly poorer fit to the data. 
The velocity amplitudes and rotations periods of individual 
spots were found, however, to be nearly identical (to within 
the 1-a error limits) for both values of u. This confirms that 
the results are insensitive to any error in the adopted value 
of the limb-darkening coefficient. 

As was mentioned in Section ^ above, the equivalent 
width of a spot bump is related to the area of the spot. 
To a very good approximation, the equivalent width of an 
absorption line in the specific intensity profile from any part 
of the stellar disc is the same as the equivalent width of 
the rotationally-broadened line profile considered as a whole. 
The ratio of the equivalent width of a small starspot bump 
at disc centre to that of the entire line profile is therefore 
just the ratio of the fiux "missing" in the spot {Fspot — 
n Ii){Rspot / dY) to the total flux from the remaining limb- 
darkened photosphere, FpUot = nIo{R*/d)^ (l~u/S) — Fspot- 
Here Jo is the specific intensity of the photosphere at disc 
centre, Rspot and R^, are the radii of the spot and the star 
respectively, and d is the distance to the star. Hence 

Wspot ^ ^ f Rspot \ ^ 

Wtot ~ l-u/3 \ R, ) 

for a small, isolated spot viewed at the centre of the stellar 
disc. 

In reality the spots are neither isolated nor completely 



Starspot rotation rates 7 



AB Dor 1996 thresh = 108 phi = 0.5 to 1.0 



50 100 

Velocity (km/sec) 

Figure 5. Rotation period P versus rotational velocity amplitude 
K for candidate starspots in phase range 0.50 to 1.00. Note that 
the spot properties fall into two main groups. The main cluster 
forms a sloping sequence delineating the stellar surface differential 
rotation. The group at top right are produced by aliasing between 
closely-spaced spots observed near the beginning and end of the 
run. The solid curve shows the fitted differential rotation law. 



at similar latitudes, with periods of order 0.5135 day with a 
phase difference of 0.02 to 0.03, observed 4 to 6 days apart. 
Aliases will occur when the data set is phased at periods 
which cause a genuine spot signature seen near the start of 
the run to appear superimposed on that of a neighbouring 
spot at similar latitude observed near the end of the run. 
Aliases will thus occur at periods given by 



alias 



1 



1 



true 



beat 



A phase shift of 0.03 rotations over a 4-day interval corre- 
sponds to a beat period of 133 days. Hence we expect a pair 
of spots with P = 0.5135 days separated by 0.03 in phase to 
produce aliases at periods of 0.5155 and 0.5115 days. Pairs 
of spots with separations of 0.02 to 0.04 rotations, observed 
4 to 6 days apart, can thus produce aliases with periods 
anywhere between 0.515 and 0.5165 days. This corresponds 
closely to the period range in which the aliases are seen in 
Fig. pi The weak candidate spot at ~ 38 km s~^ and 
P ~ 0.5134 days appears to be another type of alias, arising 
from crosstalk between the blue-shifted and red-shifted ends 
of two distinct starspot trails from low-latitude features at 
phases 0.58 and 0.67. 



dark. At the red wavelengths observed, the continuum sur- 
face brightness of the spots is expected to be about 0.25 to 
0.3 times the photospheric value. The bump amplitude for 
a spot of given area is thus expected to be diminished by 
this amount relative to the bump amplitude for a completely 
dark spot. However, the fractional coverage of spots on AB 
Dor is probably of order 30 percent, which serves to increase 
the loss of light from an individual spot as a fraction of the 
total amount of "clean" photosphere. Overall, the effects of 
finite spot surface brightness and a high spot filling factor 
will tend to cancel each other out. 

The spots listed in Table |l| have equivalent widths rang- 
ing from W = 0.0016 to = 0.0071 km s"^ Since the 
equivalent width of the deconvolved stellar profile is approx- 
imately 4.5 km s~^ (cf. Fig. |l|), the inferred fractional areas 
{Rspot/Ri,f of individual spots range from 0.00026 to 0.0012 
of the visible hemisphere. AB Dor has a radius similar to the 
Sun, so if they are considered as single, circular spots, their 
radii range from 11300 km up to 24000 km. The smallest 
detectable spots on AB Dor are thus comparable with the 
largest observed sunspots, whose umbral radii occasionally 
exceed 10000 km. 



4.2 Aliasing between neighbouring spots 

The relationship between the velocity amplitudes K and ro- 
tation periods P of the candidate spots is plotted in Fig. ^ 
Although many of the spots lie on a sloping sequence de- 
lineating the stellar surface differential rotation pattern, a 
second cluster of points with high values of K and anoma- 
lously long periods is also present. 

The upper cluster of points appears to result from alias- 
ing between closely-spaced pairs of low-latitude spots. The 
top panel of Fig. ^ shows numerous low-latitude spots, with 
a typical separation in phase of between 0.02 and 0.04 rota- 
tions. The candidates with K > 80 km s~^ and P > 0.515 
day are invariably sandwiched between a pair of candidates 



4.3 Fitting a differential rotation law 

JJonati &: (Jollier (Jameron (lyyYj found that a differential 
rotation law of the form 



fib sin 9 



(2) 



with Qeq = 12.2434 rad day"^ and fit = 0.0564 rad day"^ 
gave a good description of the latitude dependence of the 
surface rotation rate of the magnetic features and dark 
starspots on AB Dor in 1995 December. This corresponds to 
an equatorial rotation period .^equator ~ 

27r/!^e9 = 0.5132 

days, and an equator-pole lap time -P^gat ~ 27r/57i, = 111 
days. 

Combining eqs. (0) and ^ gives 

K = vsmicos6{l — sin^ 9), (3) 

which can be solved iteratively to obtain the stellar latitude 
9 at which a spot will have a projected rotational velocity 
amplitude K, given the values of v sin i, Qeq and The ro- 
tation frequency Q of the spot (and hence its period) follows 
from eq. (H). 

The equatorial rotation speed is known accurately from 
many previous Doppler imaging studies to be v sin i = 91± 1 
km s~^. We can therefore use eqs. and (^) to determine 
the values of fleq and Qt that give the best (least-squares) 
fit to the data. 

Fig. 1^ shows contours of bls a function of -Pgquator 
and -Pjjgat- Note that, since the scatter in the period residu- 
als is much greater than the precision of the measurements, 
an additional source of variance must be added to the errors 
associated with the periods listed in Table ^. As discussed 
in Section 4.4 below, the additional variance appears to be 



intrinsic to the star, and is equivalent to an additional RMS 
uncertainty of ±0.00017 d in the period determination for 
each spot. Adding this error in quadrature to the measure- 
ment errors reduced the for the fit such that its minimum 
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Figure 6. Contours of x ^ function of equatorial rota- 
tion period and equator-pole lap time. Both periods are ex- 
pressed in days. The inner contour is the locus of values = 
^min ^^"^'^ extremities in Pgquator and P^eat gi™ the 

one-dimensional Icr error bars on the two periods. The outer con- 
tours, at = . -I- 2.3 and = y2 _,_ g 2 contain 68.3% 

' ^ '^mm ^ '^mm ' 

and 95.4% of the joint probability respectively. 



value was of order the number of degrees of freedom (in this 
case, 18 data points minus 2 fitted parameters.) 

This fit yields Pgquator = 0.51319 ± 0.00007 days, and 
^beat ~ ^ ^ days. The parameters f'equator ^beat 
are strongly correlated. The differential rotat ion parameters 



for the 1995 Decemb er dataset, obtained by Donati & Col 



lier Caneron (1997), lie within the inner contour, close to 
the major axis of the error ellipse, indicating excellent agree- 
ment between the two methods. A re-analysis of th_e 1996 
Decem ber da ta using the sheared-image method of 



Donati 



et al. (|000j gives Pgquator = 0.51332 ± 0.00012 days, and 



^beat ~ ^ days, again in good agreement. 



4.4 Residuals and spot size 

The scatter about the fit to the data shown in Fig. ^ is con- 
siderably greater than the sizes of the errors on the data 
points. The measured periods of some of the best-observed 
spots deviate from the mean relation by as much as 0.005 
day, giving an accumulated error of order an hour over the 
6-day span of the data. This is too great to be ascribed to 
any plausible source of systematic error, such as irregularly- 
shaped spot groups spanning a range of latitudes. We are 
forced to conclude that the scatter in rotation rates for spots 
at any given latitude must be intrinsic to the spots them- 
selves. 

Fig. may provide a clue to the physical origin of this 
scatter. Here we plot the residuals in rotational velocity 
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Figure 7. Residual rotational velocities (expressed in m s~^) 
from the fitted differential rotation law, plotted as a function 
of spot equivalent width W. Note the increased scatter in the 
residuals for the smaller spots. 



measured relative to the mean fit, as a function of spot size. 
The scatter about the mean relation appears to be consid- 
erably greater for the small spots than for the larger spots. 
This size-dependent scatter in the rotation rates of individ- 
ual spots suggests that the smaller spots may be more sus- 
ceptible to random buffeting by su pergranular flows (see, e.g. 
Howard, Gilman & Gilman 1984). There does not, however, 



appear to be any significant difference between the average 
rotation rate of the large spots (i.e. those having W > 0.003 
km s~^), and that of the smaller spots. 



5 CONCLUSIONS 

The results of this investigation suggest that spectral 
subtraction, least-squares deconvolution and matched-filter 
analysis provide a powerful diagnostic method for identify- 
ing large numbers of individual starspot signatures in trailed 
spectrograms of AB Doradus and, by implication, other sim- 
ilar stars. The matched-filter analysis allows the positions 
and rotation periods of individual spots to be determined 
to high precision, and so provides a powerful alternative to 
methods involving regularised image reconstruction. 

We have successfully measured the rotation periods of 
individual spots at latitudes ranging from the stellar equator 
up to a latitude of 55°. The differential rotation law fitted 
to the larger spots sh ows excellent agreement with the m ean 
relation obtained by Donati & Collier Cameron (1997). The 



rate of shear in the surface rotation speed as a function of 
latitude appears to be very similar to the Sun, with the 
equator pulling one full turn ahead of the pole every 110 
days or so. 

The rotation rates of the smaller spots at the low and 
intermediate latitudes studied here appear to have a larger 
intrinsic scatter than the large spots. Further observations 
of similar or better precision would be desirable to verify 
whether this size dependence is present. If confirmed, this 
result could give a first glimpse of the turbulent sub-surface 
velocity field in a late-type star other than the Sun. 
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